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Abstract: The process parameters of enzymatic hydrolysis and molecular weight distribution 
of glutamine (Gin) peptides from soybean meal were investigated. The Protamex® 
hydrolysis pH of 6.10, temperature of 56.78 °C, enzyme to substrate ratio (E/S) of 1.90 and 
hydrolysis time of 10.72 h were found to be the optimal conditions by response surface 
methodology (RSM) for a maximal degree of hydrolysis (DH) value of 16.63% and Gin 
peptides content at 5.95 mmol/L. The soybean meal was hydro lyzed by a combination of 
Protamex® and trypsinase so that DH and Gin peptides would reach 22.02% and 6.05 mmol/mL, 
respectively. The results of size exclusion chromatography indicated that the relative 
proportion of the molecular weight <1000 Da fraction increased with DH values from 
6.76%, 11.13%, 17.89% to 22.02%, most notably the 132-500 Da fractions of hydrolysates 
were 42.14%, 46.57%, 58.44% and 69.65%. High DH values did not lead to high Gin 
peptides content of the hydro lysate but to the low molecular weight Gin peptides. 

Keywords: glutamine peptide; soybean meal; enzymatic hydrolysis; response surface 
methodology; molecular weight distribution 
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1. Introduction 

Soybean is one of the most important oilseeds in the world. The industrial processing of soybean 
yields two products, oil and a low-valued defatted soybean meal. Soybean meal protein has shown 
markedly higher values and is good source of bioactive peptides [1,2]. Glutamine is the most abundant 
amino acid in the human body, where it plays a number of important physiological roles that fuel the 
rapid proliferation of cells (fibroblasts, immune and gastrointestinal cells) [3]. The glutamic acid/glutamine 
content in wheat, corn and soybean meal is 31.90%, 20.48% and 18.82%, respectively [4-6]. Glutamine 
deficiency may compromise many important cellular protective, immunosuppressive and disease 
recovery processes that affect patient survival [7]. Prabhu reported that oral supplementation of 
glutamine or glutamic acid altered the brush border membrane in the intestine to prevent postoperative 
complications after surgical manipulation [8]. 

Recent studies suggest that plasma glutamine levels are a useful indicator of immunocompetence 
and overtraining syndrome and assert that glutamine is a potentially useful dietary supplement for 
subjects engaging in intensive athletic training [9]. However, the solubility and stability of free 
glutamine were found to be limited, and it is easily broke down into pyroglutamic acid and ammonia 
during autoclaving [10], which rules out incorporation into currently available nutritional preparations. 
Following this, the dipeptide concept with stable, highly soluble synthetic glutamine was investigated 
to overcome these drawbacks. Pep tide -bound glutamine extracted from plant protein by enzymatic 
hydrolysis has the advantage of good safety in enteral nutrition. Wheat gluten treated with protease 
was developed to produce a peptide-bound glutamine for enteral nutrition and the oligopeptide fraction 
enhanced mucosal protein content significantly more effectively in fasting rats when compared to a 
simulated free amino acid mixture. Additionally, rats showed greater resistance to methotrexate-induced 
enterocolitis when they were fed the oligopeptide fraction rather than a simulated free amino acid 
mixture [11]. Dietary supplemented glutamine and soybean peptides have important effects on the 
structure of the kidney filtration barrier and the content of total protein, albumin and 2-microglobulin 
in the urine of an over- training group [12]. Plant protein hydro lysates were mainly used as protein 
ingredients or supplements in food or as ingredients in special formulations for clinical nutrition [13]. 
Molecular weight plays an important role in the functional and biological properties of peptides [14], 
and therefore, it is essential to determine molecular weight distribution of protein hydrolysates and the 
process parameters that control it. This research aimed to optimize glutamine peptide production from 
soybean meal and to analyze the molecular weight distribution of hydrolysates, since they promise 
better utilization of soybean by-products from extracting oil. 

2. Results and Discussion 

2.1. Statistical Analysis and Model Fit of Hydrolysis Process 

The central composite design of the four independent variables (pH, temperature, E/S and reaction 
time) and the responses (DH value and Gin content) are shown in Table 1 . Multiple regression analysis 
was performed on the experimental data to obtain the predictive model using the following 
second-order polynomial Equations (1,2): 
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Ri:DH = 15.257 - 0.795 Xi + 0.118 X 2 + 0.651 X 3 + 0.379 X4 - 0.589 XiX 2 + 0.197 XX 3 - 0.048 
X1X4 - 0.227 X 2 X 3 - 0.092 X2X4 - 0.013 X3X4 - 0.313 Xi 2 - 0.359 X 2 2 - 0.200 X 3 2 - 0.139 X4 2 

R 2 :Glncontent = 5.647 - 0.075 Xi - 0.006 X 2 + 0.066 X 3 + 0.085 X4 + 0.002 XjX 2 + 0.111 X1X3 - 
O.O53X1X4 + 0.004 X 2 X 3 - 0.004 X 2 X4 + 0.079 X3X4 - 0.062 Xj 2 - 0.019 X 2 2 + 0.001 X 3 2 + 0.026 X4 2 

Table 1. Four-factor central composite design and values of responses. Independent 
variables: Temperature (Xi), pH (X 2 ), E/S (X 3 ) and reaction time (X 4 ). 
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Significance of the coefficients to DH was determined using the F-test and p-value for the influence 
of pH, temperature, E/S and reaction time (Table 2). The ANOVA of the quadratic regression model 
shows high significance for the model (p < 0.0001) and it is therefore suitable for monitoring 
optimization. The effect of pH, E/S and reaction time exerted highly significant effect on DH values 
(p < 0.0001), and hydrolysis temperature was also significant (p < 0.05). Moreover, the quadratic term 
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of pH (Xi 2 ), temperature (X 2 2 ), E/S (X 3 2 ) and reaction time (X 4 2 ) also showed highly significant effects 
on DH values (p < 0.001). In addition, two-way interaction for parameters was significant for pH and 
temperature (p < 0.0001), pH and E/S, and temperature and E/S (p < 0.001), and not significant for 
interaction between pH and reaction time, temperature and reaction time or E/S and reaction time 
(p > 0.05). The statistical analysis for the model (Table 2) showed the lack of fit was not significant for 

2 2 

pure error (p > 0.05). The coefficient of determination R was 0.9818 and Adj R was 0.9649, i.e., the 
model fit could explain 96.49% of total variability within the range of values studied. In further 
analysis, each observed value for DH (Y\) was compared with the adequate predicted value (Yi) 
(Figure 1 A). All these results imply that the model gives a satisfactory mathematical description of the 
hydrolysis process Equation (1). 

Table 2. Significance of regression coefficients for degree of hydrolysis (DH) and 
glutamine (Gin) content. 
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Significant at p < 0.05; Significant at p < 0.01; Significant at p < 0.001; Significant at 



p < 0.0001; ns = not significant. 
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Figure 1. Predicted response vs. actual 
glutamine (Gin) content. 
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The ANOVA of the quadratic regression model gives no significance for the Gin content model 
ip > 0.05) (Table 2). These variables had no significant effect on response (p > 0.05). Given R 2 at 
0.6255 and Adj R 2 at 0.2760, the model gave a poor fit with the experimental data and a disappointing 
mathematical description for Gin yield. Figure IB shows that the actual gap between Gin response (Y2) 
values and the predicted values (Y 2 '). In general, peptides activity is affected considerably by the DH 
of the protein substrate(s), and the independent variables (pH, temperature, E/S and reaction time) have 
an effect on the DH. The high peptides activity was not only a result of extensive hydrolysis, but also 
due to some intrinsic properties of the protein substrate(s) or specificity of the enzyme(s) [15]. The Gin 
peptides content was not observed to change obviously with increasing DH, so the Gin content was not 
appropriate as an objective function. The molecular weight of peptides greatly affects their functional 
and biological properties and the low molecular weight of Gin peptide results in faster absorption inside the 
body, so the activity of peptide segments were affected by the DH value to some extent [2,14]. 

2.2. Effect ofpH, Temperature, E/S Ratio and Reaction Time on Hydrolysis 

Response surface methodology (RSM) is an effective, widely accepted method of solving 
multivariate problems to properly predict the values of the response variables in many types of 
research [16]. It generates 3D response surface and 2D contour plots to show the interrelationships 
between two tested variables and the relationship between responses and experimental levels of each 
variable. The optimal values of the selected variables were obtained by regression analysis on 
Design-Expert 7.0. Two variables within the experimental range are depicted with the third variable 
kept constant at zero. Different shapes of contour plots indicate different interactions between two 
variables and circular plots of response surfaces suggest negligible interaction between corresponding 
variables while elliptical or saddle-shaped plots highlight significant interaction between 
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corresponding variables. For saddle contour plots, the optimum values are obtained at the point of 
intersection of lines formed by joining the locus [17]. 



Figure 2. Response surface plots and contour plots for the interactive effects of variables 
on the DH value, (a) pH and temperature; (b) pH and E/S ratio; (c) pH and reaction time; 
(d) temperature and E/S ratio; (e) temperature and reaction time; (f) E/S ratio and reaction time. 
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The 3D surface and 2D contour plots (Figure 2) were drawn to illustrate the principal interactive 
effects of the independent variables on the dependent DH. Figure 2A shows the 3D surface and 
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contour plots of the effect of pH and temperature on DH. Visibly, DH increased with pH until a peak at 
about 7.0. There was an ellipse in the contour plot for DH expressed as a function of pH and 
temperature, which indicates the significance of interactions between pH and temperature 
(Figure 2A). In the E/S range of 1 .0%— 3.0%, considered as the economical usage range of the enzyme, 
the DH value increased with the E/S and peaked at about 2.5% (Figure 2B). The saddle contour plots 
in Figure 2B plainly show that interactions between pH and E/S ratio also reached their optimum 
levels. DH barely increased with reaction time from 9 to 11 h as seen in Figure 2C. Low pH favored 
the higher DH value, which dropped linearly as pH increased from 6.3 to 8.0. The circular contour 
plots of response surfaces suggest that the interaction of pH and reaction time is negligible with the 
corresponding variables (Figure 2C). Figure 2D shows that DH increased gradually as temperature 
rose in the medium until gradually peaking at about 49 °C. DH also tended to increase along with E/S. 
The elliptical contour plots indicate the significance of the interaction between temperature and E/S 
(Figure 2D). As shown in Figure 2E, DH increased until the temperature reached an optimal peak at 
about 50 °C. In terms of reaction time, DH slowly increased as hydrolysis was prolonged from 9 to 1 1 h. 
Circular contour plots, which signal non-significance for the interaction between temperature and 
reaction time, were clearly shown. Further, Figure 2F shows that DH increased along with E/S from 1% 
up to 2.8% and peaked at about 2.8%. Likewise, DH rose slightly as reaction time increased from 9 to 
1 1 h. The plot in Figure 2F for DH as a function of E/S and reaction time was circular, thereby 
illustrating that this interaction was not significant. 

2.3. Optimization and Validation of the Experimental Design 

Optimal conditions of the variables were pH 6.10, 56.78 °C, E/S 1.90 and a reaction time of 10.72 h 
as determined on Design-Expert, under which the predicted DH value was 16.78%. Applying these 
conditions, the experiment scored 16.63% while Gin content was 5.95 mmol/L. These two percentages 
are close enough to validate the model. The optimum parameters of trypsinase previously reported by 
our laboratory were: pH 7.83, 50 °C, E/S 4.5 and hydrolysis time 3.7 h such that the DH value and Gin 
could reach 12.79%) with glutamine content at 5.92 mmol/mL [18]. The defatted soybean meal was 
hydrolyzed by trypsinase and Protamex® at their respective optimal parameters such that DH and Gin 
reached 22.02%> and 6.05 mmol/mL, respectively, so the soy meal was well hydrolyzed. The 
diversification of hydrolysis degree and molecular weight distribution were investigated during the 
process of defatted soybean meal hydrolysis. 

2.4. Effect of the Diversification of Hydrolysis Degree on the Molecular Weight Distribution 

The molecular weight distribution profiles of the defatted soybean meal hydrolysates obtained with 
trypsinase and Protamex® are presented in Figure 3. The molecular weight distribution of defatted 
soybean meal hydrolysates ranged from 100 to 10,000 Da and most soybean meal peptides were 
smaller than 1000 Da (Figure 3). According to their molecular weight ranges: >10,000, 10,000-5000, 
5000-2000, 2000-1000, 1000-500, 500-132 and <132, the seven absorption peaks in the 
chromatogram were marked out. Comparing the relative areas, there could be some differences in 
interpretation of the areas due to the lack of the molar extinction coefficients required to do an accurate 
quantitation. The surface area of the >1000 Da group decreased while that of <1000 Da groups 



Int. J. Mol. Sci. 2012, 13 



7490 



increased using four DH values ranging from 6.76%, 11.13%, 17.89% and 22.02%, respectively. The 
relative proportions of molecular weight <1000 Da fraction of the hydrolysates were 71.96%, 74.55%, 
86.52% and 92.37% at DH values of 6.76%, 11.13%, 17.89% and 22.02%, respectively. The relative 
proportions for 132-500 Da fraction of hydrolysates were 42.14%, 46.57%, 58.44% and 69.65% with 
DH values of 6.76%, 11.13%, 17.89% and 22.02%, respectively. The results indicate that defatted 
soybean meal hydrolysates contained many short peptides, greatly due to trypsinase and Protamex®. 
The functional and biological properties of protein hydrolysates correlated with the degree of 
hydrolysis and molecular weight [14,19]. 

Figure 3. Effect of DH on the molecular weight distribution of hydrolysates. 
(A): DH = 6.76%; (B): DH = 11.13%; (C): DH = 17.98%; (D): DH = 22.02%. 
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3. Experimental Section 

3.1. Materials 

Defatted soybean meal was supplied by Xuchang Bangdi Protein Industry Co., Ltd. (Henan, China). 
Protamex" (40,000 U/g) and trypsinase (48,000 U/g) commercial enzymes were kindly donated by 
Novozyme (Beijing, China). L-Gly-Gln standards were purchased from Sigma Chemical Co. (St. Louis, 
MO, USA). All other reagents used were of analytical grade. [Bis(trifluoroacetoxy)iodo]benzene was 
purchased from Shanghai Yurlic Chemical S&T Co.Ltd (Shanghai, China). Isothiocyanic acid phenyl 
ester was purchased from Acros Organics (Pittsburgh, USA). Acetonitrile (chromatographically) was 
purchased from Tianjing Kemiou chemreagen Co.Ltd (Tianjing, China). Ethylic acid, natrium aceticum, 
pyridine, triethylamine, trifluoroacetic acid (TFA) and formaldehyde were of analytical-grade. 

3.2. Enzymatic Hydrolysis 

The defatted soybean meal concentration was fixed at 12% (w/v) and hydro lyzed with Protamex® 
according to the hydrolysis conditions defined by the experimental design. Hydrolysis pH was 
maintained at the desired value by continuous addition of 0.1 N NaOH and 0.1 N HC1. After the 
required digestion time, the reaction was stopped by heating the solution to 80 °C for 20 min to 
deactivate the enzyme. DH and Gin content were determined from the supernatant after centrifugation 
at 4000 x g for 15 min [20]. 

3.3. Determination of the Degree of Hydrolysis 

DH was the percentage of cleaved peptide bonds (h) out of total number of such bonds in the 
substrate (h tot ) and was calculated from the amount of base consumed as given below 
(Equation (3)) [2,21]: 

h x 100 B xN n x 100 

DH% = — = -5 — (3) 

htot a x M p x htot K ' 

where B is base consumption in mL; Nb is normality of the base; a is average degree of dissociation of 
the a-NH 2 groups; M v is mass of protein (N x 6.25) in g; h is the hydrolysis equivalents in meqv/g 
protein and h tot is total number of peptide bonds in the protein substrate (7.75 meqv/g soy protein). 

3.4. Quantitative Analysis ofGlutamine in Soybean Meal Hydroly sates 

Conversion of glutamine residues into diaminobutyric acid (DAB A) was performed by adding 100 |xL 
bis-l,l-trifluoroacetoxy-iodobenzene (BTI) into acetonitrile (10 mg/mL), and 25 |xL aqueous pyridine 
(50 |j,mol/mL) was added into 100 |xL hydrolysis protein solution and incubated at 50 °C for 20 h. This 
prehydrolytic reaction was generated for glutamine residues with BTI and converted glutamine to the 
corresponding diaminobutyric acid (DABA) and was then dried immediately by vacuum. 

Acid hydrolysis was performed by adding 200 |xL of 6 M HC1 into 200 jxL DABA sample. The 
obtained solutions were purged with nitrogen gas for 2 min, then hydrolysis was carried out at 
1 10 °C for 23 h. The acid digested samples were freeze dried. 
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HPLC analysis was performed on a Waters 2695 Alliance System (Waters Corporation, Millipore, 
Milford, MA, USA), UV detector and a Waters symmetry column CI 8 oven set to 40 °C, flow rate 
0.3 mL/min and injection volume 5 |xL. The freeze dried acid digested samples was dissolved in mobile 
phase and loaded on a 2.1 x 50 mm CI 8 column (Waters Corporation, Millipore, Milford, MA, USA), then 
separated and eluted using an aqueous sodium acetate-acetonitrile gradient and detected at 254 nm [22]. 

3.5. Size Exclusion Chromatography 

Molecular weight distribution of peptides in the different hydrolysates was determined by gel 
permeation chromatography (TSKgel 2000 SWXL 300 mm x 7.8 mm) as described by Dong [23]. The 
sample was dissolved in mobile phase and filtered by micropore membrane and the injection volume 
was 20 |xL. The liquid chromatographic system consisted of a Waters 600 automated gradient 
controller pump (Waters Corporation, Millipore, Milford, MA, USA) and a 2487 UV detector at 
220 nm. The mobile phase isocratic elution consisted of 0.1% (v/v) TFA and acetonitrile (55:45). The 
flow rate was 0.5 mL/min. System control and data processing was performed using Empower GPC 
software (Version 2.0, Waters Corporation, Millipore: Milford, MA, USA). A molecular weight 
calibration curve was prepared from the elution time of the peak volume using five standards: Cytochrome 
C (12,500 Da), aprotinin (6500 Da), bacitracin (1450 Da), Gly-Gly-Tyr-Arg (451 Da) and Gly-Gly-Gly 
(189 Da). A relationship between the retention time and the log of the molecular mass of peptides used 
as standards was established. In the different hydrolysates, peptides were sorted in 5 fractions covering 
the ranges of 0-500 Da (fraction V), 500-1000 Da (fraction IV), 1000-2000 Da (fraction III), 
2000-5000 Da (fraction II) and above 5000 Da (fraction I). The relative areas of each fraction were 
given in percent of the total area. 

3.6. Experiment Design 

Temperature (Xi), pH (X 2 ), E/S (X 3 ) and reaction time (X4) were chosen as independent variables 
and optimized using a central composite rotatable design (CCRD) [13]. Four key independent variables 
at five levels were retained and individually coded as -2, -1, 0, +1, +2 (Table 3). The parameters and 
their ranges were chosen on the basis of the preliminary experimentation data not shown. All 
experiments were done in quadruplicate, results were averaged and are presented as Yi and Y 2 for DH 
and Gin, respectively. The behaviour of the system is explained by the following Equation (4): 

3 3 2 3 

y = Po + Z p£Xi + Z p " Xt2 + Z Z ^ iiXiXj (4) 

i=l i=l i = 1 j = i + 1 

where y is the dependent variable (DH value and Gin content); (3 0 is a constant; [3„ [3„ and p,y are 
coefficients estimated by the model; and x,-, and Xj are levels of the independent variables. They 
represent the linear, quadratic and cross product effects of the Xi, X 2 , X 3 and X 4 factors respectively. 
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Table 3. Coded and uncoded setting of the process parameters. 



Process Parameter 


Code 






Level 








-2 


-1 


0 


+1 


+2 


pH 




6.0 


6.5 


7.0 


7.5 


8.0 


Temperature/°C 


x 2 


40 


45 


50 


55 


60 


E/S (%) 


x 3 


1 


1.5 


2 


2.5 


3.0 


Reaction time (h) 


x 4 


8 


9 


10 


11 


12 



3. 7. Statistical Analysis 

Design Expert 7.0 statistical software (Stat-Ease Inc., Minneapolis, MO, USA) was used to analyze 
the experiment design. One-way variance analysis (ANOVA) (p < 0.05) was conducted for the 
response values obtained by the RSM model. 

4. Conclusions 

The effects of four independent variables on the production of glutamine peptides from soybean 
meal were determined using response surface methodology as a predictive tool. The parameters of the 
DH model were estimated by multiple linear regression, obtaining a good fit with the experimental 
data, given that the Adj R 2 was 0.9649. With a combination of Protamex® and trypsinase, DH values 
and Gin peptides content were 22.02% and 6.05 mmol/mL, respectively. At DH values of 6.76%, 
11.13%, 17.89% and 22.02%, the relative proportions of molecular weight <1000 Da fraction of 
hydrolysates were 71.96%, 74.55%, 86.52% and 92.37%; more notably the 132-500 Da fractions were 
42.14%, 46.57%, 58.44% and 69.65%. The hydrolysates of defatted soybean meal contained many 
short peptides and were greatly degraded by trypsinase and Protamex®. 
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